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A Weak Gravitational Lensing Analysis of Abell 2390 
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a^ : 

■ ABSTRACT 

pL| I We report on the detection of dark matter in the cluster Abell 2390 using 

O I the weak gravitational distortion of background galaxies. We find that the 

^ I cluster light and total mass distributions are quite similar over an angular 

^ I scale of ~ 7' (l/i~^Mpc). The cluster galaxy and mass distributions are 

1/^ I centered on the cluster cD galaxy and exhibit elliptical isocontours in the 

^ I central ~ 2' (280/i^^kpc). The major axis of the ellipticity is aligned with the 

CN I direction defined by the cluster cD and a "straight arc" located ^ 38" to the 

^ ' northwest. We determined the radial mass-to-light profile for this cluster and 

^ ' found a constant value of (320 ±90)/;. Mq/Lqv, which is consistent with other 

' published determinations. We also compared our weak lensing azimuthally 

I ' averaged radial mass profile with a spherical mass model proposed by the 

^ I CNOC group on the basis of their detailed dynamical study of the cluster. 

^ ■ We find good agreement between the two profiles, although there are weak 



indications that the CNOC density profile may be falling more steeply for 
e>3' (420/i~ikpc). 
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1. Introduction 



Understanding the nature of the mass distribution in clusters of galaxies is a key 
cosmological issue. Typically, three independent methods to probe the cluster mass 
distribution have been applied: 1) applications of the virial theorem to the cluster galaxy 
velocity distribution, 2) analyses of the X-ray emission of the diffuse, intracluster gas 
and 3) mass determinations based on the gravitational lens distortion of background 
objects by the foreground cluster. The latter technique, the analyses of both the strong 
as well as the weak gravitationally induced distortions in the images of faint background 
galaxies, is particularly attractive as it does not require any assumptions about the 
geometry or the dynamical of the clusters, with the weak distortions being especially 
suited for mapping the cluster mass distribution out to large radii. 



Over the past few years, there has been considerable activity to develop a theoretical 
framework for determining the mass distribution in the cluster based on a measured set 
of gravitationally distorted galaxy shapes (eg. [Tyson et al. 1990| ; [Kaiser fc Squires 1994| ; 
Schneider 1994[ ; [Kaiser et al. 1995| ; [Schneider fc Seitz 1995[ ; [Seitz, C, fc Schneider 1995[ ; 
Seitz, S., fc Schneider 199"5| , [Squires fc Kaiser 1995[ ). Several groups have developed 
algorithms and software designed to measure and analyze very small distortions expected 
in the cluster outskirts ([Kaiser, Squires fc Broadhurst 1995| ; [Bonnet fc Mellier 1995| ; 
Fischer fc Tyson 1996|) . These techniques have been applied to the central ~ l/i'^Mpc of 



several galaxy clusters (eg. Tyson et al. 1990 ; Fahlman et al. 1994 ; Bonnet et al. 1994 
Tyson fc Fischer 1995[ ; [Squires et al. 199"5[ ). 



For a few clusters, the results of the weak lensing analyses has compared with the 
results of X-ray analyses as well as detailed spectroscopic studies. There is no apriori 
reason to expect an agreement between the different analyses and in fact, such joint 
analyses can potentially enable to us learn much more about the relative distributions 
and the dynamical states of the gas, galaxies and dark matter in the clusters than is 
currently known. Indeed, in the central regions of some clusters (eg. A2218 and A1689), 
the mass inferred from the presence of (strong) lensing features exceeds that derived 
from the X-ray data under the standard assumption of pressure supported hydrostatic 
equilibrium by a factor of ~ 2-3 ( [Miralda-Escude fc Babul 1995[) . A similar discrepancy 
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also arises in galaxy cluster MS 1224, where the weak lensing inferred mass is ~ 2 times 
larger than the mass determined from a virial analysis of the cluster galaxies ( [Fahlman 
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I. 1994| ; |Carlberg et at. 1994| ). Also, in a joint lensing/X-ray analysis of A2218 ( [Squires| 



et al. 1995| ) , the lensing mass profile was found to be systematically a factor of ~ 2 larger 



than the X-ray mass profile, although the two are formally in agreement within the 95% 
confidence level. A recent analysis based on ASCA data confirmed this discrepancy 
( [Loewenstein 1996 ). Conversely, a similar analysis of the cluster A2163 yielded good 



agreement between the weak lensing and X-ray mass estimates ([Squires et al. 1996|) 



In this paper, we discuss the dark matter and galaxy distributions in the cluster of 
galaxies Abell 2390. Abell 2390 is an Abell richness class 1 cluster (|Abell et al. 1989|) 
at redshift z = 0.2279 with a measured line of sight velocity dispersion a = 1093 km/s 
( Parlberg et al. 1995| ). Abell 2390 has an X-ray luminosity of L^; = 9.77 x 10"^^ ergs/s 
(H = 75 km/s/Mpc) in the 2-6 keV energy range ( [Kowalski et al. 1984| ) and 



= 4.67 X 10^^ ergs/s in the 0.7-3.5 keV range QUlmer et al. 1986|) . A2390 has been 



the focus of several studies (eg. [Pello et al. 1991] , [Kassiola et al. 1992| , [Narasimha fc 



Chitre 1993| ) seeking to account for the observed lensing features such as the "straight 



arc" located approximately 38" away from the central galaxy ( [ir'ello et al. 199l| ). 

In this analysis, the total mass distribution in A2390 is determined using the weak 
gravitational distortion of faint background images induced by the cluster. The data 
reduction and analysis procedure used here is very similar to that we employed in our 
analysis of A2218 ( pquires et al. 1995|) and A2163 (^quires et al. 1996| ) and we refer to 



those papers for further detail. Here we will concentrate mainly on the results. We also 
compare our results for the morphology of the mass distribution and the estimate for the 
total mass with published results from galaxy spectroscopy and X-ray analyses. 



2. Data Acquisition and Analysis 

The cluster was observed at the 3.6m Canada-France-Hawai'i telescope on the 
nights of 1994 June 6-9. The detector used was the 2048 x 2048 Loral 3 CCD at prime 
focus with a pixel size of 0'.'207. For the cluster Abell 2390, we obtained 5x15 minute 
exposures in the I-band and 2 x 30 minute exposures in V with seeing conditions of 
FWHM=0'.'6 and 0'.'8 in I and V respectively. We observed the central 7' of the cluster, 
covering a square of side ~ l/i~^Mpc at the cluster redshift. 
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A detailed discussion of the observations and data reduction are described in ^quires 
et al. 1995| . Very briefly, each bias subtracted image was divided by a median twilight 
flat. The residual rms scatter in the sky background on the individual images was 23.9 
magnitudes per square arcsecond in I, and 25.3 magnitudes per square arcsecond in V. 
The data was calibrated against photometric standards in the globular clusters M92 
and NGC 4147 (unpublished photometry from Davis 1990; see also Stetson and Harris 
1988 and Odewahn et al. 1992) and Landolt (1992) standards in SAllO. Color terms 
were found to be unnecessary in the transformation and the I and V zero points were 
determined with a formal error of less than 0.005 mag. 

We detected objects with a significance of 4cr over the local sky background in 
each frame and selected objects that were detected on least two of the images in each 
waveband independently. We formed master catalogues of 1828 objects in I and 1255 
in V which are complete to ^ 23.5 and ~ 24.0 in I and V respectively. We measured 
the position, shape and brightness parameters for all of the objects using our standard 
procedure ( [Kaiser, Squires fc Broadhurst 1995| ). 



3. Cluster Light Distribution 

Combining the V- and I-band observations, we determined the color of objects 
detected in both wavebands. We identified a red sequence of objects with mean color of 
^ — / ~ 1.7 (uncorrected for reddening) at the bright end which we classified as cluster 
galaxies. To extract this color selected cluster galaxy sample, we fitted a linear model to 
the color sequence and selected bright objects (/ < 21) with color within 0.2 magnitudes 
of the mean. 

To obtain a qualitative description of the cluster light distribution, we calculated the 
galaxy light and galaxy surface number density. We employed the color selected sample 
and display as contour plots and 3D surface plots the cluster galaxy light distribution 
and surface number density in Figure ^. The distributions have been smoothed with 
a Gaussian filter with scale 0^5 (this scale was chosen match the smoothing used in 
the lensing analysis). The galaxy number and light distributions are very similar. The 
dominant structure in both distributions is centered on the giant elliptical galaxy and 
exhibits some ellipticity in the central regions, with the major axis oriented in the 
direction towards the the straight arc located ~ 38" to the northwest of the central 
galaxy. There is a small secondary peak approximately 4' to the northwest of the cD 
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galaxy that appears in both the number and luminosity weighted plots. The height of 
the secondary peak is ~ 20% of the amplitude of primary concentration in the smoothed 
light map. 

We estimated the total light content in the cluster in two ways. First, we formed 
a lower bound on the cluster light by only including the light associated with galaxies 
having the mean cluster galaxy color. This excludes cluster members with color outside 
the main cluster color sequence as well as galaxies with magnitudes beyond the faint-end 
limit. Secondly, we calculated the light in the cluster by adding the light associated 
with all the galaxies detected in the image. Since this sample includes both cluster 
as well as field galaxies, the resulting measure of the light provides an upper limit to 
the light content of the cluster, down to the limiting magnitude of the catalogue. We 
calculated the zero redshift V-band luminosity by Ly = iooMMvq-v+dm+k+Av) j^^^ 
where My© = 4.83 is the total solar V magnitude. We applied a K-correction of 
K = 0.5 as suggested by Coleman, Wu and Weedman (1980). Since A2390 is at redshift 
z = 0.2279, its distance modulus is DM = 39.28 (for = 1 and h = 1) and we used an 
extinction correction of Ay = 0.25. 

The cumulated light profiles are shown in Figure |1|. The solid line shows the 
results using all the galaxies in the images, while the dashed line shows the light profile 
computed from the color selected sample. The slope and amplitude of the light of the 
full galaxy sample exceeds that of the cluster selected subsample and for 9 > 100", the 
former is a factor of > 1.5 greater. The true cluster light profile is likely to lie in between 
the two estimates. Regardless, in computing quantities such as the cluster mass-to-light 
ratio, we adopt a conservative approach and use the higher of the two profiles. 



4. Lensing Analysis 

We identified galaxies as objects having half-light radius greater than 1.2 times the 
mean stellar half-light radius and applied a magnitude cut of 21 < / < 24. The bright 
magnitude limit was chosen, somewhat arbitrarily, to exclude cluster and foreground 
galaxies, but the results are very insensitive to the actual threshold level. With these 
cuts, we obtained 922 galaxies for this analysis. We present here the I-band data for 
the lensing analysis as the seeing conditions were significantly better than the V-band 
observations (although the V-band data was used as a consistency check). To convert 
the measured galaxy shapes into estimates of the cluster surface density, we corrected for 
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Fig. 1. — The cumulated V-magnitude and Ly luminosities from the central eUiptical 
galaxy. The solid lines come from adding the light from all galaxy candidate objects in 
the field. The dashed line is for the color selected sequence. 
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the effects of anisotropy in the point spread function (psf) due to wind-shake, guiding 
errors, etc., and the smoothing due to seeing using the procedure described in Kaiser, 
Squires & Broadhurst (1995). 

We used the seeing and psf-corrected galaxy shape measurements to determine the 
surface mass density in the cluster and display the results in Figure |^. The surface 
density reconstruction was determined using the maximum probability extension to the 
original Kaiser & Squires (1993) algorithm ( [Squires fc Kaiser 1995| ). We employed 15 



wavemodes and a regularization parameter of a = 0.05 although we note that the results 
are quite insensitive to variations in a by a factors of ~ 5. 

The peak of the mass concentration coincides with the central giant elliptical galaxy 
and is significant at the ^ 6a level. For 6 > 1', the mass isocontours are elongated in 
the southeast-northwest direction, in alignment with the direction defined by the central 
galaxy and the straight arc. There is an extension towards the north at a distance of 
~ 2' of the central galaxy which lies in the direction of the overdensity that appeared 
in the cluster number and light distributions. The peak itself (if it exists in the mass 
distribution) is not resolved. This is not unreasonable given the number of galaxies used 
in this reconstruction: assuming that mass scales directly with light, then we would only 
expect to detect the secondary mass fluctuation at the ~ la level. Thus, with the surface 
density of galaxies available here, we can not map with high statistical significance the 
mass associated with small fluctuations in the cluster galaxy distribution. 

To determine the total mass in the cluster, we employed the statistic (Kaiser et al. 
1995) 

((^1,^2) = 2{i-el/el)-' r dHe){^t) (4-i) 

which measures the mean dimensionless surface density interior to ^i, relative to 
the mean in an annulus 6i < 6 < 62- To convert ( into a physical quantity, we 
estimated T^cru = (4vrGD//3)^^ where Di is the angular diameter distance to the lens and 
P = max(0, (1 — wi/wg))- In an Einstein de-Sitter universe with Q = 1, the comoving 
distance w is defined as w = 1 — 1/Vl + z. We estimated P by an extrapolation of the 
faint redshift surveys ([Lilly 1993| ; [Tresse et al. 1993|) to fainter magnitudes and find, for 



A2390, Scrit = (6.5 ± 0.6) x IO^^/iMq/Mpc^. 

Using the estimate of Scrit? we convert the ( estimates into a mass. In Figure |^ 
we display the radial mass profile of Abell 2390 and the radial mass-to-light ratio. We 
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Fig. 2. — The bottom two rows show the luminosity and number weighted 
distributions of the color selected cluster galaxy sample. The top row is the 
reconstructed mass distribution. The distributions have been smoothed with a 
Gaussian filter with scale 0^5. In the contour plots, North is to the right; East 
is up. The box size is 7' = 966/i~^kpc. This figure can be obtained at 
ftp : / / magic bean, berkeley.edu /pub/ squires / a2390/ massandlight .ps.gz 



also plot the mass interior to the position of the straight arc in a bi-modal mass model 
proposed by Pierre et al. (1995) as well as display the predictions from two simple models 
for the cluster mass profile. The solid line shows the prediction for the mass profile from 
a singular isothermal sphere with the observed velocity dispersion ( |Carlberg et at. 1995| ), 
correcting for the mass contained in the control annulus at each radial bin. This mass 
profile does not appear to agree with the weak lensing mass profile, predicting little mass 
for 9 > 150", although uncertainties in the weak lensing estimates are sufficiently large 
that the singular model cannot be excluded with a high statistical significance. The 
dashed line corresponds to a model proposed on the basis of the CNOC analysis of the 
cluster galaxy velocity distribution ( |Carlberg et al. 1996|) under the assumption that the 
light traces the mass: p(r) = poa/r{r + a)^. We have projected the three-dimensional 
density model and corrected for the mass in the control annulus. The lensing results and 
the CNOC model are consistent within the 95% CL at all radii probed here. It is worth 
noting, however, that weak lensing profile is rising systematically more steeply than the 
CNOC profile for 6 > 150". 

We caution that in the very central regions, the lensing results should be interpreted 
with care. First, the analysis performed here is strictly only valid in the weak lensing 
regime; the presence of the straight arc at a radius of ~ 38" makes this assumption 
dubious in the innermost arcminute or so. Secondly, detailed analyses of other clusters 
seem to indicate that the lensing signal in the central regions {6 < 1') is most susceptible 
to dilution due to the contamination of our background galaxy catalogue by foreground 
and cluster galaxies. Since the results presented here were determined using only 922 
galaxies, any contamination would significantly degrade the lensing signal. Fortunately, 
the dilution effect does not appear to significant beyond a radius of ~ 1'. One way 
to address this is to select only galaxies redder than the cluster color sequence - these 
should mainly lie behind the cluster. We repeated the above analysis using such a sample 
and found consistent results, although the associated uncertainties are substantially 
larger since the number of galaxies in the sample is much smaller. 

Using the cluster light and the lensing mass profiles, we computed the cluster 
mass-to-light ratio as a function of radius. To account for the subtraction of the material 
in the control annulus in the mass estimator, we formed an analogous quantity, (l, for 
the light which is the mean light surface density interior to a given radius, relative to 
the mean in the control region. Under the assumption that mass traces the light — 
the cluster light and the lensing mass profiles are consistent with this assumption, the 
ratio of the ( values based on the observed shear and Cl forms an unbiased estimate 
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of the mass-to-light ratio in the cluster. The light estimate comes from using all the 
hght associated with the galaxies detected in the image (which, as we have discussed 
previously, includes contributions from galaxies not associated with the cluster) . We find 
a mass-to-light ratio of M/Ly = (320 ± im)hMQ/LvQ at ^ 2'.b 0.345/i-^Mpc, with 
the radial profile being consistent with a constant M/ Ly across the cluster. 



5. Discussion 



We have mapped the hght and mass distributions in the cluster of galaxies A2390 
over a scale of ~ 250" from the cluster center. The projected light distribution was 
determined using a color selected galaxy catalogue. The total radial light profile 
was calculated using all of galaxies detected, which we expect should include some 
contributions from galaxies not associated with the cluster. To map the mass distribution. 
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Fig. 3. — Left panel: the radial mass profile. The solid line is the prediction for an 
isothermal model with the measured velocity dispersion of cr = 1093 km/s and the dashed 
line is the CNOC model (Carlberg et al. 1996). Both predictions have been corrected for 
matter in the control annulus to facilitate direct comparison with the lensing data. The 
"x" is the mass estimated by Pierre et al. (1995) using a bimodal model. Right panel: 
the radial mass-to-light profile, determined by taking the ratio of the mean surface mass 
and light densities as a function of radius, relative to the mean in the control annulus. 
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we measured the shapes of faint galaxies which have been gravitationally distorted, 
correcting for systematic effects in the shape determinations such as psf anisotropy and 
seeing. 



The projected hght and mass distributions are displayed in Figure ^j. Both the light 
and mass distributions are centered on the cD galaxy. The light and mass isocontours 
are elliptical in the central region, with the major axis running from the southeast to 
northwest. There is a smaller, secondary concentration in both the cluster galaxy and 
light distributions located approximately 4' to the northwest of the cluster center. The 
mass map gives the impression that the mass distribution is extended in the direction 
of the secondary concentration seen in the cluster light distribution but there is no 
corresponding peak in the mass map. As we noted previously, this is not surprising given 
the number of galaxies used in constructing the mass map. 

Several authors (eg. Pello et al. 1991; Kassiola et al. 1992) have attempted to 
account for the lensing features seen in A2390 by suggesting that the cluster has bi-modal 
mass distribution, with a small secondary mass concentration located near the observed 
arc position. This model is further supported by the cluster X-ray emission, which shows 
elliptical isocontours centered on the cluster cD galaxy as well as a significant secondary 
peak located 40" to the northwest of the cluster center ( Pierre et at. 1953 ). The weak 



lensing mass map presented here is not inconsistent with the notion of a bi-modal 
mass distribution. As we have noted previously, the small-scale results from the weak 
lensing analysis are subject to greater uncertainties arising from the combined effects 
of nonlinearity in the lensing, contamination by cluster galaxies, and the statistical 
uncertainties due to the small number of galaxies used in the mass reconstructions. 

We calculated the mass-to-light ratio as a function of radius (Figure |^) and found 
a constant value of M / Ly = (320 ± 90)hMQ/ Lqv, although we have noted that the 
mass estimates may be biased downwards in the innermost ^ 1' due to contamination 
by faint cluster galaxies. In the course of their study of several intermediate redshift 
cluster, the CNOC group derived a mass-to-light ratio, K-corrected to redshift zero, 
of M/Lr = (331 ± 54)hMQ/ Lq^ for the cluster. We convert the Gunn r luminosity 
into V-band luminosity by applying the transformation V-r=0.16 (which is typical 
for a nearby early type galaxy) and find that the CNOC result corresponds to 
M/ Ly = (370 ± 60)hMQ/LQy, which is in excellent agreement with out determination. 

We furthermore compared the radial mass profile determined from this analysis 
with the CNOC model. We projected the CNOC 3D mass model and calculated the 
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2D radial mass profile, correcting for mass in the control annulus to facilitate direct 
comparison with the lensing results (Figure All of the mass estimates based on the 
lensing analysis agree within the 95% CL with the CNOC model (although we note that 
the CNOC model assumes spherical symmetry and there is some indications that this 
might be invalid on the scales probed here). At the very largest radius, there is a slight 
indication that the lensing inferred mass might be falling slower than the CNOC model, 
although the values are consistent within 2a. Larger scale observation of this cluster will 
prove useful to determine if this trend continues at large radii. 

As with many of the weak lensing studies performed to date, this analysis was 
confined to a relatively small field of view. This has the immediate consequence that, 
since the mass determinations made by this technique are differential, the masses quoted 
form a lower bound on the total mass at any radius. The future prospects for this 
method are exciting as the mass distributions are mapped out to large distances in the 
cluster halos. Indeed, the technology now exists for such large field lensing observations 
with the MOCAM 14' and the UH ~ 0.5° cameras at CFHT. With the types of 
observations possible with these instruments, the weak lensing technique offers a unique 
opportunity to probe many of the outstanding puzzles regarding the dark matter content 
and distribution in the universe. 
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